Jumps in current-voltage characteristics in disordered films 
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We argue that giant jumps of current at finite voltages observed in disordered films of /nO, TiN 
and YSi manifest a bistability caused by the overheating of electrons. One of the stable states 
is overheated and thus low-resistive, while the other, high-resistive state is heated much less by 
the same voltage. The bistability occurs provided that cooling of electrons is inefficient and the 
temperature dependence of the equilibrium resistance, R{T) is steep enough. We use experimental 
R{T) and assume phonon mechanism of the cooling taking into account its strong suppression 
by disorder. Our description of details of the I — V characteristics does not involve adjustable 
parameters and turns out to be in a quantitative agreement with the experiments. We propose 
experiments for more direct checks of this physical picture. 

PACS numbers: 72.20.Ht, 73.50.Fq, 73.63.-b 



Recent experiments on the I-V characteristics of disor- 
dered films of InO ^0 and TiN [3 , 4 (Similar findings 
were reported earlier for amorphous YSi films [5 ) are 
quite intriguing. Ohmic resistance R of these films at 
temperatures T < IK demonstrates insulating behavior 
fit by Arrhenius law {ks = 1) 



i?(T)=i?oexp[(A/Tn, 7 = 1, 



(1) 



with Rq and A :^ 1 ^ 5K being T-independent. Though 
this behavior is quite interesting by itself [6 , it is the I-V 
characteristics at low temperatures which turned out to 
be most unusual. When the voltage V is increased from 
F = 0, the current X first increases gradually, remaining 
rather small [highly resistive (HR) state]. At a certain 
voltage Vyll^ ^ jumps up by several orders of magnitude 
and a low resistive (LR) state arises. When the voltage 
is decreased from V > Fhl, an inverse jump between the 
LR and HR states occurs at a voltage Vlh < ^hl- These 
HR-LR switches persist in a wide range of magnetic fields 
B O^IOT, with the threshold voltages I^lh(hl) increas- 
ing with B. Authors of Refs. [11171 assumed that the HR 
state is a new collective state - "superinsulator" . Their 
calculation was vigorously disputed in Ref. [8 . 

In this Letter we show that these /- V characteristics 
can be explained without using new concepts. Our phe- 
nomenological approach predicts exactly this behavior 
when T-dependence of the resistance is steep as in Eq. 
and the electron-phonon (e-ph) thermalization is ineffi- 
cient at low T. We assume that (similar assumptions for 
a system without the switches were made in Ref. [9 ) 

1. The electron-electron (e-e) interaction is strong 
enough for electrons being mutually thermalized, i.e. 
one can introduce their temperature Tei although the 
system is driven out of equilibrium by a finite voltage; 

2. The e-ph interaction is weak, so that electrons can be 
out of equilibrium with phonons (or any other ther- 
mal bath) of temperature Tph, i.e. T^i > Tph; 



3. i?(T)-dependence at a finite voltage is the same as in 
the ohmic regime, but Tei(V) is substituted for T. 

These assumptions are sufficient to explain qualita- 
tively the experimental I — V characteristics, most care- 
fully investigated in the accompanying Letter ^ . 

Phenomenological analysis - Temperature Tq\ acquired 
by the electron due to a fixed external voltage V is de- 
termined by the balance between the Joule heating and 
the cooling by the phonon bath (electron-electron col- 
lisions conserve the energy and do not affect the heat 
balance): Before specifying the model for the electron- 
phonon (e-ph) coupling we analyze the cooling using 
model-independent arguments. The heat balance equa- 
tion for a sample of volume V and the electron density of 
states u (assumed to be smooth at the Fermi level) is 
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where £ = tt'^uVT'^/G. The r.h.s of Eq. Q obeys the 
following requirements: (i) it vanishes in the thermal 
equilibrium; (ii) it provides the entropy growth; (iii) at 
Tel ^ Tph, it ceases to depend on Tph as under such con- 
dition the radiation of the hot phonons by hot electrons 
dominates. Factor T/h provides the proper dimensional- 
ity for the relaxation rate. The second factor in r.h.s. of 



Eq. (2b) describes a suppression of the e-ph coupling at 
low T. As the phonons are gapless, the suppression is a 
power law. The energy scale encodes the strength of 
e-ph coupling. We will specify /3 and later. 

For a steep enough R{T), there is a region on the V- 
Tph plane (for Tph < T^l and V< < V < V>) where 
Eqs. ([2| have two stable solutions: a "hot" LH state with 
relatively high Tei and a "cold" HR state with Tqi ^ Tph- 
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Such a situation was discussed in Ref. [10] in connection 
with the many-body locahzation [TP. 

To analyze Eqs. ([T]) - ([2| we use dimensionless variables 



el, ph 



Tel, ph 



v_ 



and rewrite Eq. (2a) in the form 

F(tei,tph)=«; F{x,y)= [e^/^" (x^ 



1/2 



(3) 



(4) 



The new characteristic voltage, Vq emerging in Eq. (|3]), 
is natural to define through the electric field scale Vq/L^ 
independent of the system size: 



Vo 
L 



6aoh 
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(5) 



The first factor in the r.h.s. of Eq. ([5| is uniquely deter- 
mined as a combination of local characteristics of an iso- 
lated system, ctq and A, with a correct dimensionality. 
The problem of nonlinear dissipative transport is mean- 
ingless without coupling of the electrons with a thermal 
bath, i.e. when ^ oo. This is refiected by the second 
factor in Eq. (|5|. 

At the bistability boundaries t^^'^, the derivatives with 
respect to te\ of both sides of Eq. Q are equal, i.e: 



(6) 



This equation has two positive solutions if tph < ^ph(/^)- 
Here, the critical phonon temperature, = TpJ^/A is 



(7) 



The two solutions tl[^\tph) of Eq. see Fig.jl] define 
the region of voltages v'^ < v < where two states are 
locally stable, = F (^e^'^ph)- The "hot" branch t^^ 
turns out to only slightly depend on tphi 



iel(*ph)<iel(iph)<te'l(0) 



As the result, '^^(tph) is almost independent of 
(3 
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(8a) 



(8b) 



Contrarily, ioi [3 ^ 1^ the "cold" branch may be well 
approximated (except narrow vicinity of tp^^) by 



For the upper bistability boundary, we find 



(8c) 



(8d) 
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FIG. 1: Dependence of the dimensionless electronic temper- 
ature (a) and voltage (b) bistability boundaries on the bath 
temperature tph- Here t^i and are the boundaries for the 
HR ("cold") state and t^i and for the LR ("hot") state. 
Both plots are for /5 = 6, 7=1. 



At each voltage in the bistability interval < v < 
one of the states is metastable. As it is usual for the first 
order phase transition the voltages, where the switches 
between HR and LR states happen (I^l for HR^LR and 
Flh for LR^HR switches), are determined by kinetics 
of the decay of metastable states. Theoretical analysis of 
this decay and evaluation of Fhl,lh is beyond the scope 
of this paper. Here we can predict only their bounds 



(9) 



The difference between the local instability and the 
metastable state decay transition can be ascertained from 
the slope dX/dV near the transition, see Fig. [2] 

Finally, it is important to emphasize that Eqs. (|3| - 
([5| imply a non-trivial scaling of the bistability bounds. 



y<.>=A'^/2/<'>(Tph/A). 



(10) 



which should apply to Fhl, lh provided that the switches 
are close to (a)-type, see Fig. [2] 

The I — V characteristics is most conveniently visual- 
ized by calculating the non-linear conductance G 



J = G(F,Tph)F; G(Tph,F) = 



exp [-U\{vy 

Rq 



(11) 
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FIG. 2: The schematic I — V characteristics for the switches 
occuring through the local instability [type (a)] and through 
the decay of metastable state [type (b)]. 
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see Eq. ([T]), where tei{v) is found from Eq. Q. The I-V 
curves obtained by substituting the numerical solution of 
Eq. (|4| into Eq. ^ are plotted on Fig. [sj The HR and 
LR branches are connected by unstable branches shown 
by dotted lines. The HR^LR ("cold" to "hot" state) and 
LR-^HR ( "hot" to "cold" state) switches are depicted by 
vertical dashed lines. 

If Tph is not too close to the critical temperature ([7|), 



we can use the fact tei — tph <C tph, see Eq. (8c), to obtain 



an analytic description of G(y,Tph) for the HR state: 



eln [Gi?(Tph)] 
Ci?(Tph) ' 



1 < Ci^(Tph) < e. (12) 



Note, that the dependence is universal, i.e. it holds for 
arbitrary values of exponents 7. 

Microscopic input - To quantify the developed phe- 
nomenology we rely on conventional theory of normal 
disordered metals neither seeking a microscopic explana- 
tion for, e.g., the Arrhenius law ([T]) nor involving physics 
of the insulator-superconducting transition. The current 
jumps occur in the insulating regime T < A. However, 
the thermal balance equation is valid also for T > A and 
the e-ph coupling should be a continuous function. Thus, 
our strategy is to use the theory of the e-ph interaction 
in dirty metals [121 US] for T > A and extrapolate to the 
lower temperatures. The cooling rate is determined by 
the material mass density, p, and the transverse sound 
velocity Cg. Disorder is known to suppress the e-ph cou- 
pling if the wave length of a thermal phonon exceeds the 
electron elastic mean free path, hc^/T ^ ^, (as it does 
for films of Refs. [H l El H [5] where Cs - 3 x 10^ cm/s, 
T < 1 and ^ < lOnm). The result is [13] 



e-ph(^) 



a' 



h^pcl 




(13) 



where riei is the conduction electrons density, and = 
(37r^nei)^/^ is the Fermi momentum [Ml US]- Comparing 
Eqs. (13) and (2b) we obtain 



e 



1/3 



(14) 



Substituting Eq. (14) into Eq. ([5|, and using, at T > A, 



the Drude formula ctq = {e /h){kF£){nei/k^), we find 



eVo 
L 



■ a- 



1/2- 



(15) 



Remarkably, the elastic mean free-path i - the only 
quantity, which becomes meaningless in the insulator, 
does not enter electric field scale Vq/L. It makes us be- 
lieve that Eq. (15) derived for the metallic regime T > A 



may be extended to the insulator at T < A. 

Qualitative comparison with experiment- Several con- 
clusions of our phenomenological consideration can be 
compared directly with the experimental results: 



InX 




FIG. 3: Universal I-V characteristics for /3 = 6, 7 = 1 (going 
from the top) tph = 1.15 tp^, and tph = tph, tph = 0.90 tph and 
tph = 0.75 tph- Instable regions on the S-shaped curves are 
shown by dotted lines which join the upper and lower solid 
lines at the stability boundaries of the LR ("hot") and HR 
("cold") state, respectively. The vertical dashed lines with 
arrows show possible positions for LR^HR current jumps. 
The units of X and V are arbitrary, but the values of the 
current switches in the log-scale are universal and in a good 
agreement with the experimental values. 



1. The I-V curves in [H [2l |3l H [5] look like those on 



-0 



Fig|3]for both T > T^^ and T < T^^. At T 
the mfiection point becomes increasingly pronounced. 
At T below T^l, hysteretic HR^ LR and LR^ HR 
current jumps occur. These jumps indeed reach sev- 
eral orders of magnitude; 

At a rather large interval of voltages the dependence 
of logX on V looks like linear; 

Vlh only slightly changes with Tph, while Fhl sub- 
stantially increases when phonon temperature is re- 



duced. This is what Eqs. (8b), (8d) predict. 



4. Observed Tph-dependence of Fhl is still weaker than 
the exponential dependence (8d). This is consistent, 
however, with bound ([9| and the type (b) HR^LR 
transition (see discussion before Eq. ([9|, and Fig.[2|. 
Another possible reason for a relatively slow depen- 
dence of Vhl on Tph is a deviation from the Arrhenius 
law of Eq. 

5. The assumption that Fhl < and Vlh ~ im- 
plies the shapes of the I — V characteristics to be of 
type (a) close to LR^ HR switch and of type (b) in 
the vicinity of HR^ LR switch (see Figj2|. This is 
exactly what was observed. 

Further analysis of data of Ref. [2j uncovers an intrigu- 
ing discrepancy: the observed ratio of the nonlinear con- 
ductances G{Viil)/G{V 0), see Eq. (11), is noticeably 
bigger than the bounds of Eq. ( p!2| ) and Fig. [sj valid for 
any single scale dependence R{T). This discrepancy can 
be resolved within our approach only by involving an ex- 
tra temperature scale T* < T^^ in addition to A. While 
the origin of this scale is yet to be understood, the many- 
body localization [TT] is a possibility. 
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Quantitative discussion. The magnitude of the LR-HR 
jump is almost independent of other experimental param- 
eters. Since locations of the jumps between the upper and 
lower stability boundaries are ill-defined, we can only es- 
timate the order of magnitude. For Tph = 0.75TpJ^ the 
predicted 4-5 orders in magnitude of the current jump 
(Fig. [3| agrees reasonably with experiment. Experimen- 
tal cooling is well fitted with /3 = 6 in Eq. ([2|, in agree- 
ment with the above description of the electron-phonon 
mechanism. The critical bath temperature can be esti- 
mated from the temperature scale A ^ 1.2K: although 
the Arrhenius law does not give a very good fit, this scale 
can be rather reliably extracted from the data[2^. Then 
for 7 = 1 and /3 = 6, we find T^^ ^ O.IA ^ O.IK in 
agreement with the experimental values [2]. 

More quanitative comparison is hindered by the strong 
sensitivity of the critical temeperature and switching 
voltages to the choice of 7, see Eqs. ^ - Eq. ([8|. For 
a consistent quantitative comparison with experimental 
observation, one should fit the experimentally observed 
equilibrium resistance R{T) into the left-hand-side of the 
heat balance equation Eq. ( [2a| ) and solve the resulting 
equation numerically. This is done in the accompanying 
experimental Letter [2]. 

In conclusion, there is a number of strong evidences in 
favor of the electron overheating being the main cause of 
the giant HR^LR current jumps observed in [H El [3l IH 
[5]. Direct detection of the electron overheating through 
e.g., noise measurements would be an unambiguous proof 
of this. As the cooling is a rather slow process it looks 
plausible to perform the time resolved-studies of electron 
transport, e.g. measure the current caused by a train 
of voltage pulses as a function of the pulse duty cycle. 
Such a measurement could also shed some light on the 
kinetics of the switches. Also, the scaling relation ([To| 
is very characteristic for the overheating mechanism (it 
can be verified by tuning A, e.g., by gentle annealing). 

Finally, we emphasize the importance of our interpre- 
tation of the data [H lU [3l |4l [5] in the context of the 
general theory of the electron transport. Overheating of 
the electrons is quite usual in low resistive metals [18^. 
As to insulators (resistance far in excess of h/e'^) the 
overheating was rarely [9j considered quantitatively, be- 
cause the conventional mechanisms of the low tempera- 
ture charge transport are based on phonon-assisted hop- 
ping (see, e.g., p[9^). We fully realize that our explanation 
of the current jumps contradicts the standard picture (as 
well as the Arrhenius law Eq. ([T]), 7 = 1, though). At the 
same time, we find the arguments in favor for the over- 
heating being observed in Refs. [H El [3l [H [5] to be quite 
convincing. If our explanation is confirmed, these exper- 
iments should be considered as the first (to the best of 
our knowledge) reliable evidence of the strong overheat- 
ing in the insulating state and, thus, of the existence of 
phonon-independent transport in insulators. 
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